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Nonequilibrium thermodynamic formalism is proposed to derive the flux of grainy (bubbles- 
containing) matter, emerging in nucleation and growth process. Some power- and non-power limits, 
due to the applied potential as well as owing to basic correlations in such systems, have been dis- 
cussed. Some encouragement for such a discussion comes from the fact that the nucleation and 
growth processes studied, and their kinetics, are frequently reported in literature as self-similar 
(characteristic of algebraic correlations and laws) both in basic entity (grain; bubble) size as well as 
time scales. 



I. INTRODUCTION 



Nucleation and growth (NG) phenomena appear to be ubiquitous in materials processing, nanotechnology, metal- 
lurgy and biotechnology. Without loosing generality they can verbally be taken under common denominator termed 
most frequently recrystallization and grain growth though such an overall notion implicitly assumes that some ordering 
is present in the system Certainly, if the underlying ordered physical microstructure is not present, one may quite 
generally think about bubbles-containing systems, like soap froths, foams or certain sponge-resembling assemblies 0. 

The main signatures of such phenomena as model systems are as follows 's'l : 
(i) They take mostly place under certain constraints concerning the constancy or non-constancy of the overall system 
volume; (ii) The process relies on evolution of grains constituting a specimen in such a way that the mean radius 
of the grains (bubbles) has to grow in time whereas their total number has to decrease; (iii) An observed physical 
tendency while going over sufficient number of time intervals is a drop of the grain (bubble) surface free energy in 
such a way that a low surface free energy microstructure is eventually reached. 

As expected, there are very many quite important details that can be taken into account to provide a successful 
description of NG-processes, especially of their kinetics. They have thoroughly been presented in chap. 3. They 
either emphasize a pronounced role of topology, dealing with such NG-systems as composed of convex polyhedra, or 
they offer a similar metric and combined approach, in two-dimensional systems represented mostly by MuUins-Von 
Neumann evolution formula valid for mean area of n-sided grains |^. (A Monte Carlo model for three-dimensional 
case has also been proposed recently!^.) The latter is also proposed to split into three groups: (i) NG phenomenon, 
here of the so-called normal type 0,1^ ) a deterministic process governed by a continuity equation; (ii) NG of normal 
type as standard diffusion along grain size axis; (iii) NG of normal type as a stochastic process, considered in terms 
of grain (bubble) boundary motion, with a Gaussian noisy source, exclusively, and refs. therein. 

While listed above the basic signatures of NG-processes we have right now to add an important missing point. 
Namely, that at least within a statistical-mechanical framework, they are most successfully considered when a local 
conservation law is applied. It is mostly so on the so-called mesoscopic level of description, in which the system is 
in principle treated as macroscopic but where some prevalence of quantum-oriented effects, viz correlations, appears. 
No doubt that such a description leads to the well-known Fokker-Planck-Kolmogorov (FPK) formalism, with some 
possibilities of getting useful solutions, mostly after applying suitable, physically justified, boundary conditions and 
comfortably assuming the phase space to be semi-infinite. (Some visible difficulties in providing the solutions may 
emerge when a finite phase space is taken, an assumption provoked mostly by experimental situations J^.) Moreover, 
such a mescoscopic description has proved to be useful in nanoscale, that means, in a physical scale in which emergence 
of new materials is nowadays very desirable because of possible applications in medicine, microelectronics, pharmacy 
to mention but a few. 

A crucial point, however, in proposing a FPK-type description, with suitable initial and boundary conditions (IBCs), 
is to construct the flux in the space of grain (bubble) sizes. 

In a series of papers |9|, llCt , the corresponding matter flux has been contructed in a phenomenological way, in 
which to a standard diffusion (1st Pick's law) term a drift curvature - (capillarity) dependent term has arbitrarily 
been added just for reflecting properly the essential physics staying behind it 0,|El3- 

Recently, another more fundamental way of constructing the flux has been offered. It is based on dealing with a 
NG-system as an open thermodynamic system, and on applying the principles of nonequilibrium thermodynamics in 
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the space of grain (bubbles) sizes. This immediately leads to express the entropy production by the above mentioned 
flux as well as by the corresponding chemical potential gradient. The gradient, in turn, has to be defined by the 
chemical activity (not equal to one for non-ideal systems) and by the presumed, but physically very accepted under 
concrete circumstances, physical potential that can be chosen being either of power- or of non-power form in the grain 
(bubble) size. The procedure just described can be found elsewhere [T^ ITslll^ . 

The paper's contents is distributed as follows. In Section 2, we sketch the Louat-Mulheran-Harding (LMH) model 
describing a NG-process of uniform boudnaries (or, of normal type |^). In Sections 3 and 4, respectively, we propose a 
derivation of the flux of grainy (bubbles-containing) matter based on the nonequilibrium thermodynamics formalism 
(Sec. 3) whereas in the next section (Sec. 4) we discuss some important implications that may be considered in 
the new proposed physical context near thermodynamic equilibrium, mostly. Section 5 serves for concluding remarks 
preferentially within the framework of NG polycrystallinc evolutions P, [l^ . 



II. LOUAT-MULHERAN-HARDING (LMH) MODEL: A SKETCH 

The LMH model is a bit aside or slightly forgotten statistical-mechanical model that was originally invented by 
Louat for describing evolutions of polycrystallinc NG-agglomerates in dimension one |l7j . Inspite prescribed boundary 
conditions, it was thought of to be equivalent to the standard Einstein Random Walk (RW) though realized not in 
a position space but in the space of grain (crystallite) sizes. As in the precursor model, the diffusion coefhcient was 
assumed to be a constant since an incorporation of the grain surface magnitude, proposed later by Mulheran and 
Harding in higher dimensions, is certainly not possible in the one-dimensional space. 

Although the limitations of those statistical-mechanical (ordered) agglomeration models are obvious presently (the 
specialists often argue that the LMH-model becomes partly useful while describing some evolutions in nanoscale what 
seems likely to be its biggest advantage), it is still worth facing its main ideas. 

The clue idea is that the original LMH-model [TgL IT^ is a conservative model in which a flux of grainy matter is 
supposed to be of the form 

J{v, t) = -av^-'fiv, t) - Dv^ — f{v, t), (1) 

ov 

which means that it is decomposed into two parts. Here a and D are surface tension as well as diffusion reference 
constants, respectively, and generally a ^ D holds; the independent variables v and t will be defined below, a stands 
for a dimensionality d (where d=l,2,3,...) dependent quantity, namely 

a.l-i^ (2) 

The first part of the right-hand side (r.h.s.) of Eq. iQ]) is a drift term, whereas the second part of the r.h.s. of 
Eq. stands for a diffusion term. The former is a capillarity term, proportional to the grain curvature, and 
strictly related with the well-known Laplace-Kelvin- Young law, so much explored for bubbles-containing systems 
0, 1^ Il8| in which a pressure difference appears to be a main driving force. The latter is simply the first Fick's 
law, and emphasizes a proportionality of the flux to the grain surface magnitude. (Such a term, as it stands, may 
constitute the overall flux alone: This observation has extensively been explored to describe certain agglomerations 
under negligible surface or line tension regime, i.e. explicitly without the curvature-containing drift term. The system 
becomes, however, unconstrained as far as its total volume is concerned 19j.) Both the terms incorporated in |^ are 
based on phenomenological laws, so is also the construction of the mixed convection-diffusion flux l(T|l. 
After constructing the flux, we have to plug it in the well-known continuity equation 

|/(M) = -|;J(M), (3) 

where v is the volume of a grain, f{v,t) is the distribution function of the grains at time t (being of the meaning of 
the number density), that means, f{v,t)dv is a relative number of grains of size in the volume range [v,v + dv]. 
The overall conceptual and mathematical construction is to be completed by the ICS. The suitable IC looks like 

f(v,t^O)^Foiv), (4) 

where Fo{v) is a given initial distribution of grains; in many studies 0, 0, 0, 0| Fo{v) was assumed to be of 
delta-Dirac form but some other realizations are also possible Q- ^ general conclusion is, however, that the process 
in question does not depend upon an IC, cf. 9], and refs. therein. After a sufficient time spane being overcome it 
completely forgets its initial state which is a case quite often met in many open thermodynamic systems. 
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This is not the case of the BCs, in turn. The process substantially depends upon BCs chosen. A choice that remains 
as most explored is the choice of zero Dirichlet BCs, namely 

/(«-0,t) = /(i; = oo,t) = 0. (5) 

(Notice that the LMH- model does work in a semi-infinite phase space, v S [0, cxd]; it may lead to some inconsistency, 
because volume of the individual grain can be arbitrary large and can even be distinctly larger than the overall volume 
of the system as a whole, cf. discussion in [sLIT^.l 

There is some physics staying behind the BC: The grains of zero as well as of infinite size have zero account for the 
process studied. This physical constraint is sometimes called a normality condition, and the NG-process is supposed 
to be normal when the above is true. Otherwise, one may term the agglomeration abnormal. 

For the LMH model, which is formally the case of u = a in Eq. ^ the average total volume of the whole system 
remains constant in time and takes on a finite value. This is another constraint of the LMH-agglomeration, namely, 
that it takes place under constant total volume condition. This is, by the way, not the case of other evolutions with 
negligible surface tension effect, where capillary forces do not drive the evolving assembly T^. In such a system, as is 
ascertained above, we have to get rid of the curvature dependent first term of the r.h.s. of Eq. Moreover, realize 
that we practically may discard the curvature-dependent term if we enter sufficiently mature growth stage, presumed 
that the grain radii take then on quite large values so that their curvatures become very small, cf. |^ 

What remains to be done when the problem is formulated (see the above equations) is just to solve it. It was recently 
done by means of variables' separation method not only for original LMH-system but also for its physically interesting 
variations (modifications) both in finie as well as infinite phase spaces, and can be found elsewhere ^ XO, 1.1, 16, 18]. 
Some important physical quantities of the evolving agglomerate have been calculated based on the first three statistical 
moments of the process. Then zeroth moment has always been related to the average number of grains in the system 
(it mostly suffers from an algebraic drop with time for t » 1), and the first to its overall volume which is a constant 
value. From them the mean grain radius can be obtained [Tol[llllT6| : It behaves powerly with time t in an asymptotic 
time limit, and the growth exponent is l/{d+ 1). These compose the main characteristics of the model LMH-process 
just described. (From the second and first moments a useful information about volume fluctuations can additionaly 
be taken [HEl.) 

In the subsequent section, we will be still interested in a construction procedure of the flux J{v, t), hopefully arriving 
at a form resembling that of Eq. (QJ . But another clearly non-phenomenological method will be applied therein though 
we will be not capable of avoiding phenomenology completely, see the discussion concerning the physical potential (p 
below. We will, however, be able of speculate firmly and in a physically reasonable way on the basic mechanism of the 
agglomeration in the spirit of incorporating (or, having not incorporated) the suitable entropic barrier as geometrical 
one |20| . here inevitably associated with the grain curvature. 



III. NONEQUILIBRIUM THERMODYNAMICS DERIVATION OF LMH-LIKE GRAIN GROWTH 

MODEL 

Assume as above that the continuity equation Q is the evolution equation for the agglomerating system but the 
fiux J{v, t) remains unspecified. Assume, additionaly, that there is a method of constructing it. The method is called 
mesoscopic nonequilibrium thermodyiiamics |l2l [l3 . Il4| , and starts with the Gibbs equation which represents the 
entropy virtual change, namely jl2lll4j 

5S = -^J fi{v,t)5fdv, (6) 

where / = f{v,t), T is the temperature, and fj.{v,t) is the generally time-dependent chemical potential in w-space. 
The latter is given by (/i = /i(u, t) for brevity is taken) 

^i = kBT ln{af), (7) 

where a = a(v) is an activity coefficient expressing the fact that the system is not ideal; for a = 1 it would be ideal. 
(Notice that the In-iorm of the chemical potential suggests the system were ideal but, fortunately, the condition a ^ 1 
also violates this ideality.) fc^ is the Boltzmann constant. Next, a is given in terms of a potential that for the reason 
of differentiation against the chemical potential /i let us call the physical potential, ^. a reads now [T^IT^ 

a = e*/^-«^. (8) 
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Taking the temporal derivative in Eq. © and performing a partial integration (assuming, however, that J = J{v,t) 
vanishes at the ends of the phase space), one arrives at the entropy production, a, in u-space 



1 dfi 



from which we infer the expression for the current 

J 



1 w 



We have assumed that the process is local in -y. In view of Eqs ((Tn|l . Q and l|SJl we obtain 



J = -—L{v) 



df M 
ov ov 



(9) 



(10) 



(11) 



Let us define the mobility b{v) as b{v) — yjL{v) — -^^v", where D and a may have their role played in accordance 
with what has been stated in the preceding section, where they have been specified; this will somehow suggest us the 
mechanism for the entropic barrier mentioned above (a trial toward self-similarity will also be considered), and its 
connection with the physical potential </>. The obtained flux J reads finally 



J 



d 

ov 



ov 



(12) 



The expression derived, Eq. (jl2|l . is quite general (inspite of the power-law form assumed in the Onsager coefficient 
L(y)). The main physics to play now is to propose a valuable expresion for the physical potential — 4>{v). This 
is, however, a matter of proposing a basic agglomeration mechanism for the system to evolve. It will be thoroughly 
discussed in the next section by considering power- as well as non-power forms in v. 



IV. PHYSICAL IMPLICATIONS: POWER AND NON-POWER LAWS APPLIED 

In ^3 it was shown that diffusion as well as mobility functions are very related to each other. They are also related 
with the Onsager coefficient L{v) via the density function /, or also via an inclusion of the temperature T as is in the 
case of mobility, and does not apply for the diffusion coefficient. 

The diffusion coefficient (in its full form denoted by D{v,t)) is generally assumed to be defined by the Green-Kubo 
correlation formula, in which, as integrands, the random parts, J''(w,i), of the flux (|12|l of grainy matter, in two 
suitable time instants, are involved. It reads 

Div,t)^^L{v)^Dv", (13) 



that means that algebraic correlations are assumed in the correlator based on J^(v,t) |l3lll4l|. 

Since it has been argued |^ (see, literature therein), and experimentally proved that the NG-processes, namely 
recrystallization as well as grain growth are self-similar in many respects and when based on many versatile measures 
of self-similarity both in space (grain size) and time, it would be useful to test this assertion. 

Our first proposal is to assume, inspite of the algebraic (self-similar-in-time) correlations of the flux J'''{v,t) also 
that the physical potential is of a power form 

where (j)o is a constant (it may depend upon temperature, and as suggests the construction of the Lennard- Jones (LJ) 
potential, it should be nonpositive at least, or negative even; note that the attractive part of the LJ-potcntial has to 
be taken with minus sign, and we need the Van der Waals-like attraction for preserving agglomeration), and e is a 
positive exponent. 

We have now to put -^^jip = ii^^o Eq. H12f) . and the result is as follows 

J=-Dv''^f~av^f, (15) 
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where x = « — e — 1 and cr = — (/loe^-^- This way, we get the grainy flux or strictly speaking (|12() . in a very 
satisfactory form. A direct comparison with Q from Sec. 2 suggests to let e ^ (x ^ a — 1) which yields — > const 
(very weak forces as the Van der Waals typically are!). Thus, the flux H15|l is really in the form we await very much 
except that the very fact that the surface tension reference constant a unfortunately goes to zero, see above. This 
way, the drift term is completely washed out, and we get no drift effect on the agglomeration. Thus, an entropic 
barrier due to curvature (bear in mind that oc where i? is a grain radius, and when a is taken from Eq. 
cannot be proposed based on the power-like form of the physical potential cf). This is, qualitatively speaking, the case 
considered in 19] and mentioned before. 

In contrast to the above, there is a possibility of including such a geometrical barrier when a non-power (here, 
logarithmic) form of the potential is presumed. The way of doing that is worth exploring since no one can exaggerate 
the role of curvature in physics, especially in agglomeration processes 18, 22] . The first correction in curvature, e.g. 
due to the surface tension and named usually twice the mean curvature, is realized to be the main driving force (the 
whole thermodynamic context assumes the Gibbs- Thomson or capillary length to be a characteristic length of surface 
tension shrinking action) of some oversimplified gra in g rowth mechanism |2^ whereas the so-called second correction, 
frequently related with the Tolman length |l8l l22l l24l ] and named the Gaussian curvature, shows up certain elastic 
(bending, and insistance to bending termed often rigidity) properties of an agglomerate in a time instant chosen, and 
when going over its evolution stages, in which such characteristics can certainly change. 

Thus, the proposed non-power form of the physical potential is of a logarithmic form, namely 

(/> = (/>o ln{v/vo), (16) 

where Vo > stands for the initial grain volume. Amazingly, exactly this form suits very well to what we want to get. 
The entire rationale we try to reveal thoroughly throughout the present section is that we wish to have in the drift 
term, present in either Eq. or Eq. (|12|l . a factor of v"^^ which is explicitely a clear signature of the curvature, 
being a reciprocal of the grain radius R. Thus, after simple algebra -^(f) — mutatis mutandis, we are able to 
recover the desired form of the drift term in flux and this way both the fluxes, 1^ and H12I) take firmly on the 
same form, presumed that a is that dimensionality- dependent fractional number taken from (|2Jl. To have a really 
entropic barrier an additional presumption on explicit dependence of the prefactor (po upon the thermal energy ksT 
is needed, as applied for star polymer complex solutions, for example [23, carrying however less about that their 
logarithmic part of the potential is distance-dependent, whereas ours explicitely not, cf. Eq. 1)16(1 . (We may have 
some excuse here by realizing that v cx R'^, where d is the space dimension, and R can likely better mimic a distance, 
or implicitly a position measure than v does.) 

Another issue that can be addressed is of more fundamental nature. It concerns our analogy between the RW in a 
position space and the walk along the grain size axis, see beginning of Sec. 2. There are some experimental evidences 
pQ] as well as theoretical predictions 0,^3 (look at refs. therein) that the walk in grain-size space would belong to a 
broad class of geometric Brownian motion, the probability distribution of which is a logarithmic Gaussian Thus, 
instead of v in the mathematical form of Gaussian distribution a ln{v) must appear. Therefore such a potential seems 
legitimate too. 

Accepting that the key problem, because of the RW-analogy in both the spaces mentioned, does not totally rely 
on whether we have ?; or i? in potential function argument, i.e. a distance measure in (|16|) . we may at least for a 
two-dimensional case, invoke another experimental study with a logarithmic potential. The study is on a confined 
mesoscopic system in which the so-called Wigner islands appear. These are charged balls interacting via an elec- 
trostatic potential |28|. They represent the vortices in type-I superconducting systems. It turns out that just the 
logarithmic interaction potential causes maximum compaction of the charged balls of milimeter (macroscopic) size, 
typically moving and agglomerating on a conductor plane '28^ . This is then probably no surprise that the Zn-potential 
but not that power-like (|14|1 gives more chances for the system to evolve in a constant total volume regime as 
the above mentioned LMH-system does 0, Hi [2] but another one does not _19j. This is, as is stated in Sec. 2, the 
case when the governing mechanism based on curvature driving dominates, and in which a surface (line) tension effect 
cannot be avoided. 



V. CONCLUDING REMARKS 

The concluding remarks are juxtaposed in the following way: 

1. We have applied the nonequilibrium mesoscopic thermodynamics in order to justify a phenomenological con- 
struction of the grainy matter flux given in Sec. 3 by Eq. (|12|l whereas in Sec. 2 by Eq. (Q. The fundamental 
procedure offered by nonequilibrium theromodynamics method stresses mostly the role of correlations or grainy 
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matter fluctuations in the random part of the flux, as well as emphasizes a choice of the physical potential 
applied El El El. Since we know that, inspite of many additional sub- mechanisms manifesting in the course 
of agglomeration (misorientation effects; pinning of grain boundaries; dislocation density gradients' action, seg- 
regation of molecules, and alike capillary forces are those governing the system, and we may be glad 
because of finding a "compaction" ^n-potential which recovers perfectly the basic form of the flux Such a 
proposal is to our knowledge original and new. (It is worth noticing that the compaction effect, mostly forcing 
the system to exist under constant volume circumstances, where however the overall grains' evolution scenario 
proceeds in a random manner, would quite closely resemble the close-packing effect so well pronounced in amor- 
phous as well as polycrystalline materials under evolution j3] . A signature thereof is realized to be the value of 
the growth exponent which is l/{d + 1), as is mentioned in Sec. 2, and what is characteristic of the LMH-like 
systems 0,EllEaEal-) This could be anticipated as another success of the theory that before described quite 
exhaustively homogeneous [T^ as well as heterogeneous complex nucleations jl^, se also Fig. 1 for sketch 
of the conception of random close-packing. 

2. As was stated in Sec. 2, the LMH-system always evolves under a constant volume condition, therefore a pivotal 
role of surface tension in it is absolutely no surprise (the parameter cr in or in (|12l) must clearly be nonzero) . 
Therefore also the mentioned involvement of the curvature term is needed just in order to have a visible nonzero 
effect in the drift term what ist quite important when t >> 1, that means, for a mature stage of evolution when 
the grain radii are typically big but curvatures are rather small. 

3. Applyingthe nonequilibrium thermodynamics method to a relatively simple and soluble 0, IT(1 ITll | LMH- or 
alike 0i 01 systems opens up new possibilities to understand that the geometrical-physical "object", 
which the curvature is, may very likely stand as a solely picked up entropic barrier, thou gh a task remains to 
find the form of the prefactor in H16() in exact form and being temperature-dependent [20ll25l | . 

4. The forms of the proposed physical potentials are offered neither in unique nor in a systematic ways: A choice 
of potentials known in physics is very broad 29] . But we have proposed two of them which reveals some physics 
of the NG-processes with non-negligible (LMH) as well as negligible (interacting clusters formation [Ijj) surface 
tension effects. 

5. The RW-analogy emphasized so much could be an advantage; it may, however, yield also some 'uncertaintes', 
questionmarks, as for example the exact form of the physical potential and how to go out from the position or 
distance space to the grain volume or size space. Another drawback appears to be that even the form of the 
evolution equation no matter whether with or without the curvature-dependent term, is not unique: Van 
Kampen, for example, proposes a third form, in which the diffusion function is placed in front of the Laplace 
operator on the r.h.s of Eq. Q 3QJ. Therefore, generally speaking, the local conservation law may also be 
questioned. 

6. Except the agglomeration in polycrystalline model materials very often mentioned throughout the paper, one 
can expect a certain appreciable application of the offered modeling in another phenomenon, termed dropwise 
coalescence, in which the Laplace-Kelvin- Young law (see. Sec. 2), containing the capillarity-dependent drift 
term but with another cr, modified by the ratio of the satured vapor as well as liquid densities, respectively, 
is to be taken into account to describe readily the droplets' agglomeration j^, presumably with constantly 
increasing agglomerate's total volume E3- 
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FIG. 1: Random close-packing or matter compaction effect realized in systems evolving under a constraint of constant total 
volume. The nearest (and non-intersecting) neighborhood of a shadowed middle grain chosen is composed always of d -I- 1 
nearest neighbors in a d-dimensional Euclidean space; here, a demo on the spheroidal grains is proposed for simplicity. 



